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Abstract 

The sensitivity of the lightest supersymmetric particle relic density calculation to the 
variation of the cosmological expansion rate before nucleosynthesis is discussed. We show 
that such a modification, even extremely modest and with no consequence on the cosmo- 
logical observations, can greatly enhance the calculated relic density, and therefore change 
the constraints on the SUSY parameter space drastically We illustrate this variation in 
two examples of SUSY models, and show that it is unsafe to use the lower bound of the 
WMAP limits in order to constrain super symmetry. We therefore suggest to use only the 
upper value QdmIt 2 < 0.135. 

PACS numbers: 11.30.Pb, 12.60.Jv, 95.35.+d, 14.80.Ly 



During the past decade, supersymmetry (SUSY), as one of the most promising candidates for 
new physics beyond the Standard Model, has been the focus of intensive phenomenological 
studies. Huge efforts have been carried out in order to constrain the supersymmetric param- 
eter space. Among the most powerful observables for this purpose, in addition to the direct 
searches at LEP and Tevatron, stand the WMAP limits on the relic density and B physics 
constraints. 

In this letter, we present a new analysis of the relic density constraints on the Minimal Super- 
symmetric extension of the Standard Model (MSSM), and we focus in particular on bounds 
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on two gravity mediated supersymmetry breaking scenarios, namely minimal Supergravity 
(mSUGRA) and the Non- Universal Higgs Mass framework (NUHM) in which the boundary 
conditions at high scales reduce the number of free parameters of the MSSM, allowing feasible 
phenomenological studies. 

The recent observations of the WMAP satellite [1], combined with other cosmological data, 
give evidence for the presence of a cosmological matter-like density representing about 27% 
of the total density of the Universe. The remaining 73% reveal the presence of the so-called 
dark energy. From the total matter density observed by WMAP [TJ and the baryon density 
indicated by Big-Bang nucleosynthesis (BBN) [3], including theoretical uncertainties, the dark 
matter density range at 95% C.L. can be deduced: 

0.094 < tt DM h 2 < 0.135 , (1) 

where h is the reduced Hubble constant. In the following, we also refer to the older range [3] 
which admits a larger interval: 

0.1 < n DM h 2 < 0.3 . (2) 

The lightest supersymmetric particle (LSP), provided it is stable and electrically neutral, consti- 
tutes the favorite candidate for non-baryonic dark matter. The stability requirement is fulfilled 
when R-parity is conserved, and scenarios such as mSUGRA or NUHM provide us with a LSP 
satisfying the WMAP relic density constraints [3]. 

The great accuracy of the WMAP data can therefore be used to constrain the supersymmetric 
parameters, provided the relic density is calculated precisely. The computation of the relic 
density has been realized within the standard model of cosmology [5], and implemented in 
automatic codes, such as MicrOMEGAs [6] and DarkSUSY [TJ. 

However, in the standard model of cosmology, the nature of the dark energy and the evolution 
of the Universe in the pre-BBN era remain unclear. The BBN era is the oldest period in the 
cosmological evolution when reliable constraints are derived, for temperatures of about 1 MeV. 
Successful BBN models predict that radiation was the dominant energy at that time, but no 
claim is made for much higher temperatures. In fact, in models like quintessence [8], /c-essence 
[9] or dark fluid [10], dark energy could play a role before BBN, since its density could be much 
higher at such temperatures, as was especially underlined within the quintessence model in [TTJ . 
Also, some extra-dimension theories predict negative effective energies in the Early Universe, 
which can modify the relic density [L2j . 

Therefore, the standard model of cosmology could be more complex than what we think in the 
primordial Universe, and the pre-BBN era could have experienced a slower or faster expansion. 
Such a modified expansion, even though still compatible with the BBN or the WMAP results, 
changes the LSP freeze-out time and the amount of relic density. 

To model the effects of such a modified expansion in the pre-BBN era, we add to the radi- 
ation density a new dark density, varying with temperature as 



2 



where we choose To = 1 MeV and is a constant parametrizing the density behavior. Such a 
density evolution characterizes a fluid in adiabatic expansion with a constant equation of state 
wd = Pd/pd, where Pp is the pressure of the fluid: for n D = 3 (wd = 0) the dark density 
evolves as a matter density; for n D = 4 [wr> = 1/3) as a radiation density; and for n D = 6 
(wd = 1) as the density of a real scalar field (e.g. a quintessence field) with a dominating 
kinetic term [8J. no > 6 can arise for example in extra-dimension models. We introduce the 
parameter 

_ Pd(Tq) . v 

K ° = ~n W\ ' (4) 

where p rac j is the radiation density, evolving as 

Prad( T ) = 9e&(T)^T 4 ■ (5) 

9e& * s ^ e e ff ec ^i ve number of degrees of freedom of the radiation, parametrizes the temper- 
ature at which the dark density dominates the expansion, i.e. Pd(T) > P T qa(T); the larger k d 
is, the earlier the dark density dominates. In particular, if kd = 1, the radiation and the dark 
component will be co-dominant at BBN time. Thus, imposing the radiation density to remain 
dominant at BBN time and later leads to 

n D > 4 and \k d \ < 1 . (6) 

For a usual scalar field < n^, < 6, but with a modified kinetic term can reach higher 
values. We restrict here however to rip < 8. Furthermore, if kd < 0, the extra requirement 
PD + P rac [ > should be satisfied at any time, which limits negative effective densities to rip ~ 4 
or to very low |k.d|- 

The Friedmann equation at BBN time and before reads 

H 2 = ^-(Prad + Ao) , (7) 

and the dynamics of the expansion is therefore modified, leading to a higher expansion rate 
if pd > 0, or a lower one if pa < 0. Here pn does not necessarily have to correspond to the 
density of a real component, but can be only an effective term to parametrize the modification 
of the expansion rate. 

Under the standard hypotheses, i.e. in absence of entropy production and of nonthermal 
generation of relic particles, the computation of the relic density is based on the solution of the 
evolution equation [5] 

— = -3Hn - (cr eS v)(n 2 - n| q ) , (8) 

where n is the number density of all supersymmetric particles, neq their equilibrium density, 
and (dggf ) is the thermal average of the annihilation rate of the supersymmetric particles to 
Standard Model particles. By solving this equation, the density number of supersymmetric 
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particles in the present Universe and consequently the relic density can be determined. 
We consider the ratio of the number density to the radiation entropy density, 

Y(T) = g , (9) 

with 

s(T) = h eS (T) 2 -^T 3 . (10) 

^eff * s ^ ae e ff ec ti ve number of entropic degrees of freedom of radiation. Combining Eqs. ([7j) 
and (jSJ) and defining x = m^p/T, the ratio of the LSP mass over temperature, yield 



dx V 45G x 2 

with 



9 l /2 = ^(l + J- dJ ^) ■ (12) 



v^el V 3/i eff dT 
Note that in the limit where 0, we retrieve the results of Ref. |5| 



The freeze-out temperature Tf is the temperature at which the LSP leaves the initial ther- 
mal equilibrium, i.e. T = Tf when Y(Tf) = (1 + 5)Y e q(Tf) , with 5 ~ 1.5. The relic density 
is obtained by integrating Eq. ffTTl) from x = to mj^sp/To, where T = 2.726 K is the 
temperature of the Universe today [3]: 

Q LSP h* = 2.755 x 10 8 fg|r(T ) . (13) 

To compute numerically the relic density, we use a modified version of MicrOMEGAs 2.0.7 [6] 
which includes the alteration of the expansion rate in the primordial Universe, as in Eq. (JTj). 
The SUSY mass spectrum and couplings are computed with SOFTSUSY 2.0.14 [13], and the 
b — > «7 branching ratio and isospin asymmetry are calculated with Superlso 2.0 [14] . using the 
limits of [15]. 

Let us consider first the mSUGRA parameter point (m = 600 GeV, mu2 = 730 GeV, 
Aq = —mo, tan (3 = 50, fi > 0), which is favored by WMAP, as an example. The stan- 
dard calculation of the relic density for this point leads to ^LSP = 308 GeV, ^lsp^ 2 — 0.105, 
and a freeze-out temperature of Tf = 11 GeV. In Fig. [TJ the dependence of the relic density 
and of the freeze-out temperature on no and kd is shown. The oblique black and red lines 
correspond, respectively, to the present upper limit of the WMAP constraint Q^ph 2 = 0.135, 
and to the older limit Q^gph 2 = 0.3. kd is varied in the interval [10 -10 , 1], and no in [4, 8.2]. 
First we note that the relic density can be increased by up to a factor 10 6 and the freeze-out 
temperature up to 50 GeV. This strong dependence of the relic density and freeze-out tem- 
perature on % and k d is due to the fact that for large values of n D and k , the Universe 
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Figure 1: Relic density (left) and freeze-out temperature (right) in the modified pre-BBN 
expansion model for the mSUGRA point specified in the text. The regions on the right of the 
oblique black/red lines are disfavored by the new/old constraints of (pQ) and (J2]). 



is not anymore dominated by the radiation at freeze-out temperature, but by the dark fluid. 
Avoiding the extreme values of rip and kd, and considering for example a quintessence fluid 
with ri£, = 6, we observe that even with a negligible contribution of dark fluid at BBN time, 
e.g. K£> ~ 10~ 3 , the relic density is multiplied by a factor larger than 100, and this point is then 
excluded by the WMAP limits. For lower values of no and of kd, the relic density remains 
however inside the favored interval. 

For negative effective densities the calculated relic density hardly changes. With = 4, it 
decreases by less than 1%, even with the excessive value k d = —0.5. 

We study now the effects of a quintessence density with ri£> = 6 and ko e [10~ 5 , 10~ 2 ], and of 
a density with rip = 8 and ko = 10 -5 , in the mSUGRA and NUHM parameter spaces. Note 
that this choice of parameters is completely in agreement with the known cosmological con- 
straints, and such a modification of expansion rate would be transparent on usual observables. 
Nevertheless, it would affect the constraints obtained on the SUSY parameter spaces. 

In Fig. [21 we present the constraints on the mSUGRA parameter plane (mi/2,mo), for tan/3 = 
50, A = —m and /i > 0. In the yellow region the LSP is charged, therefore this region is 
cosmologically disfavored. The green region is disfavored by the branching ratio of B — > X s 7, 
the red region by the isospin asymmetry of B — > K*j, the gray area is excluded by the collider 
constraints on the particles masses [16], and the dark blue and light blue are respectively fa- 
vored by the new and old constraints of ([1]) and (j2J). All contours are at 95% C.L. The top-left 
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plot corresponds to the relic density calculated with the usual expansion rate. From top to 
bottom and left to right k^> is increased (as well as no for the last plot), and we can notice the 
change in the relic density favored contours. 

The first important feature of this figure is that even with a modest modification of the ex- 
pansion rate, a large part of the favored zone enters into the charged LSP area. Only a very 
narrow line is still favored for = 10~ 4 — 10 -5 and no = 6, but becomes completely excluded 
for larger kd and no- Therefore, a slight alteration of the pre-BBN expansion rate can induce 
drastic modifications in the mSUGRA WMAP favored regions. 

In Fig. [3J we consider the NUHM parameter plane (//, tua), for m = 1000 GeV, m 1 / 2 = 500 
GeV, tan/3 = 35, A = and fi > with the same color definitions as in the previous figure. 
First in this plane the LSP is not charged. The WMAP contours are small strips forming a 
cross. When increasing the expansion rate alterations, the cross becomes smaller and smaller, 
and in the last plot, for rip = 8 and Kd = 10~ 5 , the region favored by WMAP is extremely 
small, yet existing. 

The analysis presented here reveals that an increase of the expansion rate before BBN, even very 
small and presently undetectable in the cosmological observations, decreases the LSP freeze- 
out temperature, and increases the relic density by a factor of up to 10 6 . Even with a tiny 
increase, the WMAP favored regions in the SUSY parameter spaces are strongly displaced, so 
that they can become excluded by other constraints. Since the pre-BBN era is relatively uncon- 
strained, and in many models beyond the cosmological standard model the expansion rate can 
be modified by the presence of a new fluid, a modified gravity, or any other reason, we consider 
that using the WMAP data to constrain SUSY should be done with caution. In particular, 
this study reveals that the relic density is increased when the expansion rate is modified by 
the presence of an extra density, but not decreased (or very slightly, in the case of a negative 
effective density). Thus, a relic density originally excluded by the lower WMAP limit could 
get increased and shifted to the permitted interval. Therefore, with such a modification of the 
cosmological expansion rate, only the upper limit of the WMAP constraints can be used safely, 
as it provides a limit on the maximum relic density. If we disregard the lower limit, we should 
not consider the region favored by WMAP anymore, but instead the region disfavored by the 
upper limit, as shown in Fig. HI In this way, two problems are avoided: first, if it turns out 
that the cosmological standard model is too simple to describe correctly the pre-BBN era and 
that the expansion rate should be modified, our conservative upper limit (in the cosmological 
standard model) would still provide valid exclusion zones; second, if dark matter is composed 
of several components - axions, dark fluids, ... - and not only of the LSP, the upper limit 
of the WMAP constraint would still be reliable, contrary to the lower limit. Therefore, as a 
conclusion, we suggest using 

tt LSF h 2 < 0.135 (14) 

in order to explore SUSY parameter spaces with the WMAP data. Also, as we can see from 
Fig. HI this WMAP constraint, together with the other usual constraints, can provide already 
good information on the favorite region of the SUSY parameter spaces. 



6 




400 800 1200 1600 2000 400 800 1200 1600 



m 1(2 [GeV] m 1(2 [GeV] 




400 800 1200 1600 2000 400 800 1200 1600 2000 

m 1(2 [GeV] m 1(2 [GeV] 



Figure 2: Constraints in the mSUGRA parameter plane for several values of and n^: rip = 0, 
Kd = (top left), n D = 6, k d = 1CT 5 (top right), n D = 6, k d = 10 -4 (middle left), n D = 6, 
k d = 1CT 3 (middle right), n D = 6, k d = 10~ 2 (bottom left), n D = 8, k d = 10~ 5 (bottom 
right). The dark and light blues correspond respectively to the regions favored by the new and 
the old dark matter constraints. The other colors refer to exclusion contours as described in 
the text. 
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Figure 3: Constraints in the NUHM parameter plane for several values of kd and no- no = 0, 
kd = (top left), no = 6, kd = 1CT 5 (top right), no = 6, kd = 10 -4 (middle left), Ud = 6, 
kd = 1CT 3 (middle right), no = 6, kd = 1CT 2 (bottom left), no — 8, Kd — 10~ 5 (bottom 
right). The color conventions are described in the text. 
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Figure 4: Constraints on the mSUGRA (left) and NUHM (right) parameter planes. Contrary 
to the precedent figures, the light blue contours stand for the regions "disfavored" by the 
upper limit of the old dark matter constraints. The dark blue region together with the light 
blue, depicts the area disfavored by the upper limit of the new WMAP results. The other 
color conventions are as in the precedent figures. The only favored zone (by all the considered 
constraints) is therefore the white area. 
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